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Abstract —In order to project electromagnetic fields between 
different meshes with respect to the conservation of energetic 
values, Galerkin projection formulations based on the energetic 
norm are developed in this communication. The proposed for¬ 
mulations are applied to an academic example. 

Index Terms —Finite element methods, Galerkin method. In¬ 
terpolation, Modeling, Projection. 

1. Introduction 

In recent years, the numerical studies of coupled problems 
are more and more investigated. These studies deal with 
the interactions between different physical phenomena, e.g. 
electromagnetic - thermal or magnetic - mechanic. According 
to the importance of the interaction, coupled problems can 
be treated with different strategies. One possibility is that the 
study domain is discretized on different meshes for different 
problems. In this case, it is necessary to communicate and 
transfer fields between different meshes. In the literature, the 
concept of Galerkin projection based on the en'or norm 
provides a very convenient tool to carry out this transfer IT), 
0, E). In comparison to the direct interpolation, Galerkin 
projection enjoys several advantages, especially in terms of 
precision. However, using this method, the conservation of 
the energy is not assured between the original and target 
meshes 11. To tackle this issue, formulations deduced from 
the minimization of the energetic norm can be used. 

In this communication, energetic approaches for Galerkin 
projection are developed in order to conserve the magnetic 
energy and electric power between different meshes. Firstly, 
the numerical models developed from the minimization of 
the energetic norm are presented for magneto-static and eddy 
current problems. Secondly, the obtained formulations are 
applied to an academic example. 

II. Energetic Galerkin projection formulations 

Given a solution of electromagnetic computation on a source 
mesh, we aim to project this result onto a target mesh which 
differs. The energetic Galerkin projection formulations are 
investigated for magneto-static problems as well as eddy 
current problems. 

A. Magneto-static problems 

In order to solve magneto-static problems, different formu¬ 
lations such as the formulations based on scalar or vectorial 
potentials can be employed. As a result, either the magnetic 


held H or the magnetic Hux density B is conformed with 
physical properties. We denote by and the helds 
obtained on the source mesh and by H, and B, the helds to 
be calculated on the target mesh. The energetic norms of the 
interpolation en'or are dehned as: 

£h= r /j|H,-H/dr (1) 

Jd 

eb= f-IB^-B/dr (2) 

JdI^ 

with p the linear magnetic permeability. 

Using Whitney elements in 3D, H, and B, are discretized 
in the edge and facet element spaces respectively such that 
H, = y and B, = y w^h,. Here (resp. wf) and /i, 
(resp. bp are the shape functions and the values of H, (resp. 
Bf) associated with the i-th edge (resp. facet). To project helds 
onto a target mesh with respect to the magnetic energy, weak 
formulations based on the minimization of the energetic norm 
are developed. In ([T]|, the energetic error norm is minimized 
when its derivatives with respect to all degrees of freedom are 
equal to zero, thus for each edge i: 


dhi 



|H, - H/ dr = 0 


X 


/j(w^ ■ H, - w^ HJdr = 0 


(3) 

(4) 


Finally, the matrix equation to solve can be written: 


[C] [h] = [F] (5) 

where Cij - I /iw- ■ wjdr, F, = I juw- ■ H^dr and [h] is the 
Jd Jd 

vector of degrees of freedom to be calculated. 

A similar demonstration can be applied to for the 
projection of B^ with respect to the magnetic energy. 


B. Eddy current problems 

For eddy current problems, either magnetic or electric 
harmonic formulations can be used in order to calculate 
harmonic helds. We obtain either Hj or conformed with 
physical properties. Thus, depending on the used formulation, 
the energetic error norm to minimize can be dehned 0: 

8h= r /J|H,-H.pt/TH- r — icurl H, - curl dr (6) 
Jd Jd, ctw 


2 


ee = r -iE,-E/fl'T+ r - 

Jd, ^ JD^t 

with cr the electrical conductivity and a> the pulsation. 

In order to project fields onto a target mesh with respect 
to the magnetic energy in D and electric power in conducting 
domain Dc, weak formulations can be obtained by the min¬ 
imization of the energetic norm. Equation (j^ gives rise to 
matrix equation: 


curl E, - curl E^ 


dr (7) 


([C] + [Ce]) [h] = [F] + [F,] (8) 


where [C] and [E] are the same matrix as in magneto-static 
problems, Ceij - J" curl w- ■ curl Wjdr and Fei = 

—curl w® ■ curl H^dr. Here [C] + [Cg] is a positive- 
, o-oj 

definite matrix. 

A similar matrix system can be obtained from (|^ in order 
to project E^. 


III. Application 



Figure 1: Geometry of the example 
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(a) Energetic Galerkin projection 
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(b) 1 } Galerkin projection 


Figure 2: Distribution of B (T) on a clipping plane, calculated 
from the projected H, field. 


The proposed projection formulations are applied in an 
academic example. This example is composed of a magnetic 
core cylinder and an excitation coil (Fig. [T]). Two meshes are 
considered: Ms (306K elements) and Mt (60K elements). In 
order to carry out projections, Ms is used as the source mesh 
for the first computation and Mt is considered as the target 
mesh for the projection. 

Firstly, a constant current is applied in the coil. The 
magneto-static problem is solved on Ms or Mt using the scalar 
potential formulation. The obtained field H, on Ms is then 
projected to Mt using 0- The corresponding magnetic flux 
density on the clipping plane (S in Fig. is presented in 
Fig. 1^. In order to illustrate the advantage of the projection 
method, the field H j is also projected using a classical l} norm 
projection (Fig. |^). In comparison to the energetic projection, 
the 1 } projection fails to provide a correct distribution of B at 
the boundary of the cylinder. Table [^presents the max values of 
fields and the magnetic energy calculated on Ms and Mt alone 
for reference, and then from Ms to Mt projections. With the 
formulation deduced from the minimization of the energetic 
norm, the magnetic energy and the max values of fields are 
close to the reference results (calculated on Mt alone). 

Secondly, a sinusoidal current is applied as excitation. In 
this case, eddy currents appear in the conductor cylinder. 
The magneto harmonic problem is solved with the magnetic 
formulation. Table [H] presents the reference results (calculated 
on Ms and Mt alone) for the max values of fields, the magnetic 
energy in the domain and the ohmic losses in the cylinder. As 
the 1 } projection fails to conserve the magnetic energy and 
the ohmic losses, only the results obtained using the energetic 
projection are shown. They are obtained using (j^. In this 
example, the magnetic energy as well as the ohmic losses are 
well conserved using the energetic Galerkin projection. 


Table I: Validation of the energetic projection approach 
(magneto-static problem) 



magnetic energy 
{mj) 

|H|inajr 

(kAjm) 

(T) 

calculated on 

Ms 

8.01 

44.1 

0.116 

calculated on 

Mt 

8.25 

49.0 

0.096 

proj. Ms- 

»Mt 

10.28 

42.9 

4.17 

energetic proj. Ms^Mt 

7.32 

45.3 

0.089 

Table II: Validation of the energetic projection approach (eddy 
current problem) 


magnetic 

ohmic losses 

^^\max 

1 JImaj: 


energy (mJ) 

(JIW) 

{kAjm) 

{kAjrrP') 

calculated on Ms 

8.05 

104 

44.9 

21.1 

calculated on Mt 

8.18 

104 

49.7 

20.2 

energetic proj. 
Ms^Mt 

7.54 

98 

46.8 

19.6 
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